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We consider dissipative dynamics of atoms in a strong standing laser wave and find a nonlinear 
dynamical effect of synchronization between center-of-mass motion and internal Rabi oscillations. 
The synchronization manifests itself in the phase space as limit cycles which may have different 
periods and riddled basins of attraction. The effect can be detected in the fluorescence spectra of 
atoms as equidistant sideband frequencies with the space between adjacent peaks to be inversely 
proportional to the value of the period of the respective limit cycle. With increasing the intensity of 
the laser field, we observe numerically cascades of bifurcations that eventually end up in settling a 
strange chaotic attractor. A broadband noise is shown to destroy a fine structure of the bifurcation 
scenario, but prominent features of period- 1 and period-3 limit cycles survive under a weak noise. 
The character of the atomic motion is analyzed with the help of the friction force whose zeroes are 
attractor or repellor points in the velocity space. We find ranges of the laser parameters where the 
atomic motion resembles a random but deterministic walking of atoms erratically jumping between 
different wells of the optical potential. Such a random walking is shown to be fractal in the sense 
that the measured characteristic of the motion, time of exit of atoms from a given space of the 
standing wave, is a complicated function that has a self-similar structure with singularities on a 
Cantor set of values of one of the control parameters. 
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I. INTRODUCTION 

Light exerts mechanical forces on neutral atoms whose 
origin can be explained as in the particle as in the wave 
pictures of the particle-wave duality. In the first one, 
light is a collection of photons carrying energy, angu- 
lar momentum, and momentum. When interacting with 
light, atoms change not only their internal electronic 
states but also their external translational states, the pro- 
cess known as a photon recoil. If one considers light as 
a wave, the mechanical force is produced by the interac- 
tion between inhomogeneous electromagnetic wave and 
the atomic dipole moment. 

In this paper we consider atoms confined in a cavity or 
in a standing wave made up of two counter-propagating 
laser waves. There is a vast literature on the various as- 
pects of the topic of mechanical action of light on atoms 
(see, for example, some reviews and books [H 0, S S ) ■ 
In the simplest case of a two-level atom interacting with 
a single mode of an one-dimensional standing wave in 
a cavity, there are, at least, three strongly coupled sub- 
systems: internal and external atomic degrees of freedom 
and field degrees of freedom. Recently, it has been shown 
theoretically and numerically that Hamiltonian dynamics 
of a single atom in a standing light wave demonstrates a 
variety of new dynamical effects such as correlations be- 
tween Rabi oscillations and atomic translational motion 
0, 0, the Doppler-Rabi resonance 0, @, Hamiltonian 
chaos 0, E3 , Levy flights El El , and atomic dynam- 
ical fractals yj fl2. . The semiclassical and quantum 
theory of these effects has been developed in the cited 



papers. 

In the present paper we study dissipative dynamics of 
two-level atoms in a standing wave pumped by a laser and 
report on a nonlinear dynamical effect of synchronization 
between internal and external atomic degrees of freedom. 
Light-induced forces cause an atom either to be trapped 
in a well of the optical potential or to fly through the 
standing wave. The mechanical atomic oscillator forces 
the Rabi one to oscillate with the same rhythm. Syn- 
chronization occurs after a transition time and manifests 
itself in the phase space as limit cycles which may have 
different periods depending on the detuning S between 
the atomic transition, uj ai and standing- wave frequen- 
cies, ujf, and on the average number of photons in the 
wave n. Synchronization and the period of the respec- 
tive limit cycle can be detected by measuring fluorescence 
spectra of atoms. Synchronized ballistic motion with a 
pcriod-m limit cycle should generate sidebands with the 
frequencies, luj ± jkfv s /m (j = 1,2,...) in the respec- 
tive spectra (where v s is a quasistationary atomic veloc- 
ity), whereas trapped atoms, oscillating under a period- 
m limit cycle, generate sidebands ujf ± jQ^/m, where 
£1^ is the respective frequency of nonlinear synchronized 
oscillations. Asynchronous chaotic motion generates a 
broaden spectrum. 

This work is organized as follows: in Sec. II the model 
of an atom interacting with a standing wave in a driven 
cavity is introduced, and semiclassical equations of mo- 
tion including spontaneous relaxation, cavity damping, 
and pumping are derived. In the strong-field limit and 
under a resonance of a laser with the standing wave, 
the equations are reduced to a five-dimensional dissipa- 
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tive nonlinear system with coupled mechanical and Bloch 
variables. In Sec. Ill we discuss briefly two simple regimes 
of atomic motion: shallow oscillations of an atom trapped 
in the optical potential and ballistic motion under an in- 
fluence of the light-induced friction force. The friction 
force is analytically calculated for thermal atoms with a 
large detuning S. In Sec. IV we report on synchronization 
of the atomic Rabi oscillations with the momentum oscil- 
lations, give respective approximated analytic solutions 
for period- 1 limit cycles, show that basins of attraction 
of different limit cycles are riddled, and compute spec- 
tra of atomic fluorescence. In Sec. V we demonstrate a 
rich scenario of bifurcations in dependence on the average 
number of photons in the standing wave. The scenario 
contains coexisting limit cycles of different periods, pe- 
riod doubling, and cascades of different bifurcations on 
the road to a strange chaotic attractor. Sec. VI contains 
numerical results demonstrating atomic dynamical frac- 
tals in the regime of random walking of an atom in a 
standing wave and a comparison between dissipative and 
Hamiltonian fractals. 



II. EQUATIONS OF MOTION 

Let us consider a two-level atom with the electronic 
transition frequency to a and mass m a in a ID-cavity 
which sustains a single cosine-like standing-wave mode 
along the axis r with the wave number kf and the fre- 
quency ujf. The cavity mode is supposed to be pumped 
by a classical laser field with the amplitude s, phase ip, 
and frequency lui. The total Hamiltonian 

H = H + Hi (1) 
consists of the standard cavity QED Hamiltonian 

2m a 2 V 2 / 

— hfl a cos(fc/f )(a t (7_ + a<7 + ) (2) 

and the laser-mode interaction Hamiltonian 

Hi = heaexpi(LUit — ip) + H. c. (3) 

The atom may emit and absorb photons from the cav- 
ity mode. Both the processes change internal atomic 
and field states. The external or translational atomic 
states are changed as well because of a photon recoil. 
The laser field supplies an energy to overcome atomic 
and mode dissipation. All the processes are presented 
in the Hamiltonians. The operators f and p with the 
commutation rule [f,p] = ih describe atomic position 
and momentum, respectively. The Pauli atomic opera- 
tors a± = \{a x ± ia y ) and <x z with the commutation 
rules [<J±,o z ] = T2(Tz, [&+,&-] = &z describe electronic 
atomic transitions. The boson operators a and with 
the commutator [a, a'] = 1 describe the cavity mode. 



The atom-field coupling is characterized by an amplitude 
of the single-photon Rabi frequency fio- 

It is convenient to choose the following combina- 
tions of operators: e = a' exp(-iuift) + acxp(iujft), 
g = i[a> exp(— iuft) — aexp(iojft)}, x = a~ cxp(ioj ft) + 
<t+ exp(— iujft), y = exp(iujft) — <r+ exp(— iujft)], a z , 
f, and p, for which the respective Heisenberg equations 
can be derived using the Hamiltonians J2J and J2J) in the 
frame rotating with the frequency u) / 
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where Tf and T a are the rate of dissipation of cavity 
photons and spontaneous emission rate, respectively, to 
be added to the respective equations phcnomcnologically. 
After averaging over an initial quantum state, which is 
supposed to be a product of the field, internal, and exter- 
nal atomic states, we get the respective c-number equa- 
tions of motion. Assuming a large number of photons in 
the cavity mode, n = (e 2 + <? 2 )/4 ^ 1, and large values 
of the atomic momentum as compared with the photon 
momentum hkf, we adopt the semiclassical approxima- 
tion. It results in the following normalized closed set of 
ODEs: 

£, = ap, 

i> = ( l /2)(gy - ea;)sin£, 

e = ycos£ — 27/e + 2£ , sin(Ar + (j)), 

g = xcos£ — 2"f f g — 2i?cos(Ar + 0), (5) 

x = Sy + zg cos £ - {l/2)^ a x, 

y = -Sx + ze cos£ - (l/2)j a y, 

z = -{gx + ey)cos£ - -f a (z + 1), 

where dot denotes differentiation with respect to the 
dimensionless time r = Hot and E = e/Hq. Here 
^ = kf<f>,z=<a z >,p=<p>/hkf, and the 
other variables are the expectation values of the respec- 
tive operators written above. A few dimensionless fre- 
quencies define the atom-field dynamics: a = fiJc 2 /m a fio 
is the photon recoil frequency, S = (ujf — cj a )/f2o the de- 
tuning between the frequencies of the cavity mode and 
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atomic electronic transition, A = (u>f— uj{)/VLq the cavity- 
laser detuning, and j a j the relaxation rates. 

In fact, Eqs. (JSJ) are the coupled Hamilton-Maxwcll- 
Bloch equations which take into account not only light- 
induced mechanical forces acting on a neutral atom but 
the back reaction of the atom on the cavity mode as well. 
In a strong field the cavity mode reaches rapidly a regime 
with a large saturation number of photons n = (E/2jf) 2 . 
Throughout the paper we will consider the case of exact 
resonance between the cavity mode and the laser, A = 0, 
and a specified phase ip = ir/4. It follows the reduced set 
of our basic equations 



£ = ap, 

p = —u sin£, 

u = 6v - (l/2)j a u, 

v = —8u + 2nzcos£ — (l/2)^f a v, 

z = — 2ucos£ — 7 a (z + 1), 



(G) 



u 



where we introduce new variables 
y/ n/2 (x + y) and v = (gx + ey)/2 



(ex - gy)/2 = 
y/n/2 (y - x) 
which arc proportional to the quadrature components of 
the atomic electric dipole moment. Namely, the system 
(jSJ) will be treated throughout the paper. Eqs. © are 
written in the form convenient to the comparison with 
equations (0) in Ref. Q which describe Hamiltonian self- 
consistent atom-field dynamics in a high-Q cavity with- 
out dissipation and a pumping field. In fact, we need no 
cavity in the strong field regime for which Eqs. JBJ) are 
valid. A standing wave may be created by two counter 
propagating laser beams in a free space. Real laser beams 
produce, of course, a three dimensional standing wave, 
and we simplify the geometry to simplify the analysis. 



III. CENTER-OF-MASS MOTION IN A 
STANDING WAVE 

A. Shallow oscillations in a well of the optical 
potential: steady-state solutions 

An atom moves in an optical potential II = J u sin£d£, 
whose form can be easily found in the steady state ap- 
proximation when the atom is assumed to move so slowly 
that the Bloch variables it, v, and z have a time to take 
steady-state values. It is the case of rather cold trapped 
atoms oscillating in wells of the optical potential that is 
valid if the normalized Dopplcr shift is much smaller than 
the detuning and the relaxation rates, \ap\ <C |<5|,7 a ,/. 
With the time derivatives of the Bloch variables in Eqs. 
(|rj|) equal to zero, it is easily to find their stationary values 

2nd cos £ 



S 2 + 2n 2 cos 2 £ + 7 2 /4' 

7 a ncosg 

S 2 + 2n 2 cos 2 £ + 7 2 /4' 



l + 2n 2 cos 2 £/(<5 2 +7 2 /4)' 



Substituting the solution for u s in the second equation 
in (JgJ), we obtain the gradient or dipole force acting on 
the atom 



P = 



nS sin 2£ 



S 2 + 2n 2 cos 2 £ + 7 2 /4 
and the optical potential 



U =-\n(^+8 2 + 2n 2 cos 2 £ 
n \ 4 



(8) 



(9) 



Well known effects follow from this solution (see, for ex- 
ample, 0, H). At exact resonance, 5 — 0, the optical 
potential is zero, and the atom moves with a constant 
velocity, equal to its initial value p = pq. At positive 
values of the detuning, 6 > 0, atoms tend to go to nodes 
of standing wave (£ = 7r/2 + irm), whereas at negative 
detuning, 5 < 0, they tend to go to antinodes (£ = irm), 
where m = 0, 1, .... We recall that our standing wave is 
of the form — cos£. 



B. Ballistic motion and the friction force 

When an atom moves with a thermal average mo- 
mentum p s through the standing wave, the Dopplcr 
shift exceeds relaxation rates, and the approximation of 
slowly varying field amplitude is no valid. In the limit 
\ap s | > 7a,/, we will seek solutions for the Bloch variables 
in the form 



u = A cos ap s r + B sin ap s r, 
v = C cos ap s r + D sin ap s r. 



(10) 



Let the atom reaches after a transition time a quasis- 
tationary momentum p s (around which its instant mo- 
mentum oscillates slightly) and does not experience Rabi 
oscillations due to a large detuning \6\ 3> 1, i. e. z ~ — 1. 
Substituting the solutions (fTU|) into Eqs. ©, it is possible 
to find after some algebra an expression for the force act- 
ing on the atom / = p = — Msin£. After time averaging 
it, we get a dissipative or friction force, F = — d\p\/dT, 
which has the form 



2nS"f a a\p s 



( la 6) 2 + [(a\ Ps \) 2 -5 2 +1 2 a /4Y 



(11) 



This approximated solution gives dispersion-like depen- 
dencies on the Doppler shift a\p s \ and the detuning 5. At 
negative detunings, the friction force decelerates atoms 
whereas at positive detunings, it accelerates them. 

In fact, the friction force F depends on a|p s | and S 
in a more complicated way. We integrate the basic Eqs. 

and plot in Fig. 1 its dependence on instant values of 
the atomic momentum \p\, averaged over a time of flight 
between two successive nodes of the standing wave, at 
S > (Fig. la) and 5 < (Fig. lb). It is plotted for 
the values of the momentum larger than a critical mo- 
mentum p cr that is a minimal value of \p\ under which 
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FIG. 1: Friction force F in units of Q.ohk[ as a function of 
the average atomic momentum |p| in units of the photon mo- 
mentum hkf. (a) Positive detuning 5 = 24 and (b) negative 
detuning S = —24. 



a ballistic motion is possible. What happens with atoms 
whose momenta at some instant of time lie in a given 
range? At blue detunings, 8 > , atoms with \p\ < \p a \ 
are trapped, where \p a | is a value of the average momen- 
tum which corresponds to the first zero of the function 
The parts of this function between its successive 
zeroes may be considered as averaged phase trajectories 
of atoms in the space (\p\, d\p\/dr). At positive values of 
the friction force, atoms move with decreasing momen- 
tum, to the left on the plot. At F < 0, they move with 
increasing momentum, i.e. to the right on the plot. Thus, 
zeroes of the function F (\p\) with positive derivatives are 
attractors, whereas zeroes with negative derivatives are 
repellors. 

As it can be seen from Fig. lb, at red detunings, 8 < 0, 
repellors and attractors interchange their places as com- 
pared with Fig. la. It immediately means that the fric- 
tion at 8 < accelerates slow atoms and decelerates fast 
ones which in a due time reach a quasistationary momen- 
tum \p s I and oscillate slightly in the velocity space around 
it (Fig. 2a). If the initial momentum is sufficiently small, 
an atom may even change the direction of motion, but a 
stabilization occurs in a due time. It is the known effect 
of velocity grouping that appears to be useful to ob- 
serve fluorescence spectra for an atomic ensemble with a 
velocity distribution that is a more feasible task from the 



experimental point of view than observing a single-atom 
fluorescence. Initially trapped atoms oscillate firstly in 
wells of the optical potential with increasing amplitudes, 
and then leave the wells reaching after a transition time 
a grouping velocity (lower trajectory in Fig. 2a). 

In difference from red detunings, where velocity group- 
ing is global, at blue detunings it is limited in the velocity 
space by a range of the atomic momenta \p a | < \p\ < \pb | , 
containing zeroes of the function F(|p|). The values of 
the momenta, \p a \ and \pb\, corresponding to the first 
and last zeroes, are proportional to the detuning <5 and 
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FIG. 2: Atomic momentum p as a function of time r in units 
of the reciprocal of the single-photon Rabi frequency fio. (a) 
red detuning, S < 0: global velocity grouping for different 
values of the initial momentum; (b) blue detunings, S > 0: 
fast atoms are accelerated (upper trajectory), slow atoms are 
decelerated (lower trajectory), and velocity grouping occurs 
only for atoms in a range \p a \ ^ \p\ ^ \pt\- 



depend on the other control parameters in an irregular 
way. Thus, at 8 > 0, slow atoms (with \p\ < \p a \ — 170 in 
Fig. la) are decelerated by the friction force and trapped 
in wells of the optical potential eventually (see the lower 
trajectory in Fig. 2b). The friction force for fast atoms 
(|p| > \Pb\ — 320) is negative, and they are acceler- 
ated by the field (the upper trajectory in Fig. 2b). The 
velocity grouping at 8 > occurs only in the range 
170 < \p\ < 320. A sample trajectory for such atoms 
is shown in Fig. 2b as a middle one. 
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IV. SYNCHRONIZATION 
A. General features 

Our basic Eqs. 10 describe two coupled oscillators, 
the external mechanical one (£,p) and the internal Rabi 
oscillator (it, v,z) whose free frequencies may differ in a 
few orders of magnitude in a strong field. The frequency 
of small oscillations of an atom in a well of the opti- 
cal potential can be estimated to be a 1 / 2 ™ 1 / 4 , whereas 
the free Rabi frequency, \/S 2 + An, is proportional to a 
square root of the average number of photons n. The 
phenomenon we report in this section is synchronization 
of internal and external degrees of freedom of an atom 
in a standing-wave field. Synchronization is a fundamen- 
tal nonlinear phenomenon to be observed for the first 
time by Ch. Huygens in 1665 with his famous mechan- 
ical pendulum clock and described theoretically by 
A. Andronov and A. Witt in 1930 [l^ in terms of limit 
cycles. In our case the effect consists in equality between 
the frequency of oscillations of the atomic internal energy 
and the frequency of variations of the atomic momentum 
as for trapped as for ballistic atoms. The mechanical 
oscillator forces the Rabi one to oscillate with the same 
rhythm. A standing wave with the spatial period 2tt 
modulates with the Doppler shift ap s momentum of an 
atom, moving with an average momentum p s . 

Synchronization manifests itself in the phase space as 
limit cycles which may have different periods. A limit cy- 
cle of period m means a periodicity occurs when the atom 
transverses 2m times the same node (being trapped in a 
well) or different nodes (a ballistic atom) of the standing 
wave. We have found limit cycles with m=l, 2, 3, 12. 
In Fig. 3a and 3b we plot projections of period-1 and 
pcriod-3 limit cycles on the Bloch plane (v, u) and on 
the plane (£,z) of the internal and external variables at 
S = 24, a = 0.01, 7 a = 0.3, and p a = 60. The size of 
a limit cycle is defined by the amplitude of the respec- 
tive synchronized oscillations, its form — by the spec- 
trum of the oscillations, and the time of circulation of a 
phase point over the cycle — by the period of oscillations. 
The larger the parameter of nonlincarity, n, the stronger 
the form of a limit cycle differs from the harmonic one. 
In Fig. 3a the dashed line is a period-1 limit cycle with 
n = 3000 and the solid one is a cycle with n = 10000. 

Synchronization of internal and external atomic de- 
grees of freedom is a common feature of regular motion 
that occurs with different initial states and in wide ranges 
of values of the control parameters. The only exception 
is the case of exact resonance, S = 0, when the optical 
potential vanishes and an atom flies through the wave 
with a constant velocity, but its internal energy oscil- 
lates with a modulated amplitude and frequency. Chaos 
destroys synchronization if one understands under this 
name a frequency locking, where two frequencies become 
rationally related. Synchronization of chaotic oscillations 
(if any) is beyond the scope of this paper. 

The time of settling limit cycles depends mainly on the 
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FIG. 3: Projections of the phase trajectories of a trapped 
atom onto two planes (u, v) and (£, z) at 8 = 24: (a) period-1 
limit cycle (dashed line at n = 3000 and solid line at n = 
10000), (b) period-3 limit cycle (n = 14846), and (c) chaotic 
strange attractor (n = 24000). The atomic position £ is shown 
in units of fcT 1 . 



action of the force / which accelerates or decelerates an 
atom during a transition time r s before it reaches a qua- 
sistationary momentum p s . This time can be estimated 
for atoms with velocities p closed to p s as follows. The 
force can be expanded in the Taylor series nearby the 
point p s 

Af = Apf(p s ) + .., (12) 

where Ap = p s -p, Af = j(p s ) - f(p), and /' = df /dp. 
On the other hand, for small Ap, Af ~ Ap/r s , and in 
the first order we get r s ~ [/'(p s )] _1 . It follows from 
(|llfl . the transition time is a reciprocal of the average 
number of photons n. 



B. Approximated analytic solutions for period-1 
cycles 

Depending on the values of the control parameters and 
the initial momentum, atoms cither are trapped in wells 
of the optical potential and oscillate with p s = 0, or they 
fly through the standing wave (p s ^ 0). The frequency 
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of synchronized motion is defined by the frequency of 
oscillations of p around p s . At S > and sufficiently 
large values of the initial momentum po, atoms can move 
with an acceleration asymptotically going to zero (Fig. la 
and Fig. 2b, the upper trajectory). 

In order to find simple analytic solutions for period- 1 
limit cycle in the case of ballistic motion, we put £ ps ap s 
and expand the variables p, u, v, and z in Fourier series 



The solutions (|16fl - (|18|) are valid under the following 



condition: n£4 <C 5 — 



72/4. 



Under the limits of 
its validicity, the solutions (| 1 fc>|) - H18|) agree well with 
numerical simulation. Since the normalized amplitude of 
small oscillations, £ rn , can be less than 1, the solutions 
(|16fl are valid up to the values of the number of photons 
of the order of a few thousands. 



3=0 



dj ' cos(jap s r) + b 



(«) 

3 



sm(jap s r), 



(13) 



with p, v, and z being analogous Fourier series with re- 
spective coefficients. Substituting these series in the ba- 
sic equations ©, we get an infinite number of recur- 
rent relations for the coefficients. A period- 1 limit cycle 
is described by the first terms in the Fourier series, if 
amplitudes of higher harmonics can be neglected. The 
approximated solution at large dctunings can be found 
(truncating the hierarchy at j = 2) to be 



P ~ Ps - (SA b /2ap s ) cos2ap s r, 

u ps — 25 Ab cos ap s r, 

v ps —jaAb cos ap s T + 2ap s Ab sin ap s T, 

z ps — 1 + Ab(l + cos2ap s T), 



where 



A,, 



(S 2 



a 2 p 2 . 



■72/4) 



(14) 



(15) 



The solutions (|14|) and (|15fl arc valid under the condition 
n <C 5 2 — a 2 p 2 + 7„/4. Numerical simulation shows that 
the obtained solutions really work well under the con- 
dition of a small number of photons as compared with 
the squared detuning. When this condition breaks down, 
higher harmonics become significant, and the simple form 
of period-1 limit cycles is distorted. 

The same procedure has been used to find analytic 
solutions for small oscillations of atoms in wells of the 
optical potential. At S > 0, the wells are situated at 
the nodes of the standing wave £ = tt/2 + irm (m = 
0, 1, ...). Assuming a harmonic motion of an atom at the 
bottom of a well, £ = 7r/2 + nm + £ m eoswjT, we find the 
approximated solutions for a period-1 limit cycles 



C. Riddling 

For given values of the recoil frequency a and the 
spontaneous rate 7 a , there are many kinds of attractors 
when the detuning S and the numbers of photons n are 
changed. For the same n and <5, different initial condi- 
tions may lead to different attractors, i. e. the system 
is sensitively dependent on initial conditions. It should 
be stressed that it may happen even under regular syn- 
chronous oscillations. We have numerically found two, 
three, and more coexisting stable limit cycles with the 
same and different periods, whose basins of attraction 
are riddled by each other. The basin is said to be riddled 
|l7j . if any point in the basin has in its arbitrary small 
vicinity points of other attractor basins. 

In Fig. 4 we illustrate riddling for three coexisting limit 
cycles of periods 1,2, and 3 in the plane of initial atomic 
momentum po an d population inversion zq. Initial con- 
ditions on 200 x 200 grid were integrated to find the re- 
spective limit cycles they destinate. Grid points going to 
the period-3 limit cycles are plotted as dark grey dots, 
those that go to the period-2 limit cycle — as light grey 
dots, and the grid points going to the period-1 limit cycle 
are left white. As one can see in Fig. 4, grey and white 
points are mixed. Since the accuracy of preparing ini- 
tial conditions is finite, one cannot predict exactly which 
limit cycle will be destinated from a given initial point in 
the phase space. To confirm the riddling, we magnified 
a small region in Fig. 4a (0.8 ^ zq ^ 1; 60 ^ po ^ 70) 
and compute the basins of attraction with the same grid 
200 x 200. The result, shown in Fig. 4b, demonstrates 
that riddled basins may occur not only on chaotic and 
strange attractors 01 but also on stable limit cycles. 



p ps (2SA w ^ m /(j() sina^r, 
u ps — 25A w ^ m cos u>£T, 
v ps -j a A w £ m coso^r + 2u£A w £ m sinter, 



(16) 



cos 



where the squared frequency of the small oscillations is 



(17) 




and 



A,„ = 



(18) 



D. Spectra of fluorescence 

The feasible way to detect synchronization in real ex- 
periments is to measure spectra of atomic fluorescence. 
Let us derive an expression for the atomic transition 
dipole moment in the laboratory frame of reference. As- 
suming the magnitude fj, of the dipole d to be real, we 
have d = p,x, where i-component of the dipole in the lab- 
oratory frame is connected to the x and y components in 
the frame, rotating with the frequency ujf, by the follow- 
ing way: 



x cos ujft + y sin ujft. 



(19) 
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FIG. 4: Riddled basins of attraction in the plane of initial 
population inversion zq and initial momentum po (in units 
hkf): □ period 1 limit cycle, □ period 2 limit cycle, □ period 
3 limit cycle, (a) 200 x 200 grids and (b) 200 x 200 grids giving 
a magnified view of a small region in (a). 



Our basic equations 10 axe written in the rotating frame 
for the dynamical variables u = \fnj2 (x + y) and v = 
sVn/2 (y — x). Therefore, the expression for the atomic 
dipolc moment in the laboratory frame is 



/' 



[(u — v) cos u>ft — (u + v) sin u tt] . (20) 



Using the solutions (|14|) for a pcriod-1 limit cycle in the 
ballistic mode, it is easily to find that the atomic dipole 
moment will oscillate with the frequencies u)f ± jkfV s , 
where v s is a quasistationary atomic velocity reached by 
atoms after a transition time, and j = 0, 1, 2, .... Dipolc 
moments of ballistic atoms with a pcriod-m limit cy- 
cle have the frequencies u>f ± jkfv s /m = ujf ± j'Slocvp s . 
Trapped atoms oscillating with a period-m limit cycle, 
have dipole moments oscillating with the frequencies 
cjf± jClg/m = ujf ± jQ,ow^/m, where lu^ is given by Eq. 



(I17fl . Fluorescence spectra both for the atoms, randomly 
walking in a standing wave and oscillating in wells in the 
regime of strange chaotic attractors, should be broad- 
ened. 
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FIG. 5: Fluorescence spectra for trapped atoms in different 
regimes of oscillations: at S = 24: (a) period- 1 limit cycle with 
n = 3000 photons, (b) period-1 limit cycle with n — 10000, (c) 
period-3 limit cycle with n = 14846, and (d) chaotic strange 
attractor with n — 24000. The frequency uj is shown in units 
of fio- 



In order to check the conclusions mentioned above we 
integrate numerically Eq. © and compute the Fourier 
spectra of the quantity (u — v) cos w ft — (u + v) sinw/t in 
different regimes of oscillations. In Fig. 5 spectra for dif- 
ferent limit cycles for the same parameters as in Fig. 3 arc 
shown. Fig. 5a demonstrate the spectrum of the quasi- 
harmonic period-1 limit cycle at n = 3000. In Fig. 5b a 
more complicated spectrum of the period-1 limit cycle at 
n = 10000 is shown. Fig. 5c demonstrate the spectrum 
of the period-3 limit cycle (n = 14846), and Fig. 5d is a 
spectrum of the chaotic strange attractor (n — 24000). 
A commom feature of all the spectra is absence of odd 
harmonics Wf ± 2jQ^/m (j = 0, 1,2, ...). For quasihar- 
monic period-1 limit cycles, it can be shown analytically 
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in the framework of the approach described in Sec. IVB, 
that odd harmonics are really absent in the oscillations 
of the dipole-moment variables. 

All the spectra refer to a single atom. A multiatom 
spectrum should be broadened because of a velocity dis- 
tribution. However, the effect of velocity grouping, dis- 
cussed in Sec. IIIB, should help to resolve the sidebands 
of the light emitted by an ensemble of atoms with limit 
cycles of low periods. 

V. BIFURCATIONS 

A. Synchronization map 

Our nonlinear dissipativc dynamical system demon- 
strates a rich variety of dynamical long-term regimes in- 
cluding periodic, quasiperiodic, and chaotic strange at- 
tractors. To illustrate this complexity we plot in Fig. 6 
a synchronization map in comparatively small ranges of 
the detuning and the number of photons with fixed val- 
ues of the other control parameters, a = 0.01, 7 a = 0.3, 
and the same initial conditions. White color corresponds 
to period- 1 limit cycles, nuances of grey color — to limit 
cycles of higher periods, and black color — to determin- 
istic chaos. As the number of photons increases, higher 
harmonics in the spectrum become more and more promi- 
nent, and the approximated solutions (|14|l - l|18|) for 
period- 1 limit cycles are no more valid. It is seen from 
the map that the limit cycles with higher periods appear 
practically for random values of the control parameters 
but inside prominent parallel bands separated by zones 
with pcriod-1 limit cycles, exclusively. The border of 
chaos in Fig. 6 is jagged, and there exist nearby small 
"islands" of chaotic motion submerged into the "sea" of 
period- 1 limit cycle motion. 

The case of chaos we identify computing the maximal 
Lyapunov exponent A which characterizes a mean veloc- 
ity of diverging of initially closed trajectories in the phase 
space. The dependence of A on the parameters 8 and n 
is shown in Fig. 7 where color modulates values of A. 
If A > the atomic oscillations are chaotic in a sense 
of sensitive dependence on initial conditions. A promi- 
nent asymmetry with respect to sign of the detuning <5 
is explained by different character of the force acting on 
atoms when 5 > and 6 < (see Sec. IIIB). 

B. Bifurcation diagram 

A general approach in studying the complexity of non- 
linear dynamical systems involves investigating qualita- 
tive changes in their dynamics as system's control param- 
eters are varied. Sudden changes in the dynamics, known 
as bifurcations, occur when one of the control parameter 
crosses a critical value. The best way to visualize it is to 
compute so-called bifurcation diagrams. The bifurcation 
diagram in Fig. 8 shows dependence of negative values of 
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FIG. 6: Synchronization map in the plane of the average num- 
bers of photons n and the detuning 8 (in units of f2o): pe- 
riod 1 limit cycle, □ period 2 limit cycle, ■ period 3 limit 
cycle, ■ limit cycles with periods from 4 to 12, and ■ chaos. 
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FIG. 7: Dependence of the maximal Lyapunov exponent A on 
the average number of photons n and the detuning S (A and 
8 are shown in units of f2o)- 



the w-component of the atomic dipole of trapped atoms 
at the moments of time when u = on the average num- 
ber of photons n at fixed initial conditions and S = 24. 

A general view in Fig. 8a demonstrates a complicated 
scenario of bifurcations on the road to chaos. In the range 
n < 12000, there exists a single stable limit cycle with 
period 1. At 12000 < n < 14000, this cycle begins to 
multiply, i.e. there coexist in the phase space a few dif- 



9 




-80 



■ ! 



(a) 



6000 



12000 



18000 



• : i§il 

^ ■ H 



24000 




17000 



18000 



in the phase space four different limit cycles of period 1 
with their own basins of attraction. A pitchfork b\ at 
n ~ 17050 means not a period doubling bifurcation but 
birth of two new limit cycles with period 1 from the old 
one. With n increasing further, real bifurcations with 
period m occur. We indicate them in Fig. 8b by the 
respective numbers. A prominent bifurcation of period 3 
occurs at n ~ 17200 when a period-3 limit cycle appears 
suddenly. Cycles of period 5 and 7 are also identified. 
Famous period doubling bifurcation are rather rare. One 
of them, corresponding to appearance a period-6 instead 
of a pcriod-3 limit cycle, is indicated in Fig. 8b as "6". 
In all the range, 16900 < n < 18000, there coexist, at 
least, three period-1 limit cycles. 
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FIG. 9: Maximal Lyapunov exponent A (in units of flo) as 
function of the average number of photons n at 5 — 24. 



FIG. 8: Attractor bifurcation diagrams at S — 24: (a) a gen- 
eral view and (b) a cascade of bifurcations. Arrows 61,3 in- 
dicate doubling of period-1 and period-3 limit cycles, respec- 
tively; 3, 5, and 7 — period-3, 5, and 7 limit cycles, respec- 
tively; 6 — period doubling bifurcation 3^6. 



ferent period-1 limit cycles. Starting with a given initial 
state, an atom after a transition time reaches, of course, 
a definite period-1 limit cycle, but small changes in initial 
state may lead to another period- 1 limit cycle. In other 
words, we know that in this range of n a period-1 limit 
cycle is settled after a transition time, but one cannot 
predict practically which one is reached. In the range 
14000 < n < 21000, a few cascades of bifurcations occur 
where limit cycles of higher periods appear. 

A more detailed view in Fig. 8b helps to see a fine 
structure of one of the bifurcation cascades. We see that 
at n < 17050 practically all the points belong to four 
parallel lines. At a given value of n, only one value of 
v exists, but under small changes in n values of v jump 
randomly between the four lines. There really coexist 



z 

(a. u.) 




FIG. 10: Strange chaotic attractor in the Bloch sphere at 
6 = 24 and n = 24000. 

It is seen from Fig. 8a that, beginning with n « 22000, 
a complicated intermittent motion takes place. To con- 
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firm it, we compute in detail the dependence of the max- 
imal Lyapunov exponent A on the average number of 
photons n at S = 24. The plot in Fig. 9 demonstrates 
an intermittent road to a "stable" chaotic mode that is 
nothing more than a strange chaotic attractor. Its view 
in the Bloch space is shown in Fig. 10 with the parameter 
values S = 24 and n = 24000. Projections of this attrac- 
tor on the planes (v,u) and (£, z) are shown in Fig. 3c. 
Its chaoticity is confirmed by computing tha maximal the 
Lyapunov exponent to be A ~ 1, and its strangeness — 
by computing the Hausdorff dimension dp w 2.7. 



C. Influence of noise 

The cascade of bifurcations in Fig. 8b is very compli- 
cated and its detailed analysis in the whole parameter 
space is practically impossible. We cannot expect, of 
course, that all the fine details of this cascade are rele- 
vant for real atoms in a standing wave. In order to know 
which limit cycles and bifurcations are expected to sur- 
vive under an inevitable noise in real experiments, we 
incorporate a stochastic force into the right side of the 
second equation in the set © and compute again the 
bifurcation diagram shown in Fig. 8b at different lev- 
els of noise. Noise destroys the fine structure but the 
prominent period- 1 and period-3 limit cycles with their 
characteristic properties do not disappear. In Fig. 11 
we demonstrate the bifurcation diagram in the range 
16900 < n < 18000 under a weak broadband noise. Esti- 
mating the frequency of nonlinear oscillations of trapped 
atoms to be in the range uj^ ~ 0.3 — 2, wc model the 
broadband noise by a linear combination of 100 har- 
monic functions with random phases and equidistantly 
distributed frequencies in the range from 0.05 to 5. 



VI. FRACTAL RANDOM WALKING OF 
ATOMS 

In Sec. Ill we have considered mechanical effects in 
regular atomic motion: acceleration, deceleration, veloc- 
ity grouping, and trapping of atoms, all of which were 
characterized by negative or zero values of the maximal 
Lyapunov exponent A. The Lyapunov map in Fig. 7 
demonstrates that there exist wide ranges of the control 
parameters, n and 5, where atomic motion is expected 
to be exponentially sensitive to small changes in initial 
conditions. In the chaotic regime there are also exist os- 
cillations in wells of the optical potential and ballistic 
flights with an average velocity p s . However, instead of 
stable limit cycles, there appear chaotic attractors (see 
Figs. 3c and 10), and the instant atomic velocity in the 
ballistic mode oscillates around p s in an irregular way. 
In Fig. 12a we plot an example of a chaotized velocity 
grouping (at S = -24, n = 24000, and p = 8) that 
should be compared with the velocity grouping in regular 
regime in Fig. 2a. In difference from the regular velocity 
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FIG. 11: Attractor bifurcation diagram at a weak noise. The 
other parameters and initial conditions are the same as in 
Fig. 8. 



grouping, the momentum p a , averaged over time of flying 
between two successive nodes of the standing wave, is not 
constant but oscillates slightly in an irregular way. With 
sufficiently large p s it does not matter, in principle, but 
in some ranges of the control parameters a cardinally new 
type of motion takes place, when atomic motion is neither 
ballistic nor oscillating in wells but resembles strongly a 
random walking. In this regime an initially trapped atom 
oscillates in a well for a while, then it leaves the well for 
a ballistic flight and can be trapped in another well, etc. 
Random walking appears when there exist such zeroes 
of the friction force F where the grouping momentum p s 
is close to the critical momentum p cr (see Sec. IIIB). In 
the chaotic regime (A > 0), p s varies in an irregular way 
being either greater or smaller than p cr . As an example, 
we demonstrate in Fig. 12b and c a random walking (at 
S = —3, n = 24000, and po — 50) with a rather small 
grouping momentum p a w 30. 

Interesting properties of randomly walking atoms will 
be demonstrated with the following numerical experi- 
ment, a scheme of which is shown in Fig. 13. Atoms, 
one by one, are placed at the point £ = with the 
same initial momentum po = 50 along the cavity axis. 
We compute the time T the atoms need to reach one of 
the detectors varying one of the control parameters, the 
average number of photons, n, or the detuning S, and 
with the other parameters and initial conditions to be the 
same. To reduce numerical efforts, wc restrict ourselves 
by ID-geometry. The analogous numerical experiments 
for Hamiltonian atomic motion in a high-Q cavity have 
been done in Ref. Q for a classical radiation field and in 
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FIG. 14: Dependence of the exit time T (in units Q, ) on the 
1000 2000 3000 4000 number of photons n with different resolutions at S = — 3 and 



la = 0.3. 



FIG. 12: Chaotic atomic motion: (a) chaotic velocity group- 
ing at 5 = —24, (b) and (c) random walking at S = —3. 
Atomic position in units kj 1 and momentum p in units /ifc/ 
as functions of time. Here n = 24000. 



atoms 




FIG. 13: Schematic diagram showing scattering of atoms at 
the standing wave when the space between detectors is equal 
to two (dashed lines) and twenty one wavelengths (bold lines). 



Refs. 0Q3 for a quantized field in the Fock and coher- 
ent states of light. Self-similar hierarchical structures in 
the dependence T(j>q) have been found in these papers. 
It seems to be more practical to compute the exit-time 
function in dependence on one of the parameters (n or 
(5) that is more easy to control than the initial atomic 
momentum. 

The exit-time function T(n) has been computed at 
S = — 3 with two detectors placed at the distance of 
two lengths of the standing wave (dashed lines in Fig. 13 
indicate places of the detectors in this case). Fig. 14 
demonstrates an intermittency of smooth intervals and 
complicated structures that cannot be resolved in prin- 



ciple no matter how large the magnification factor. The 
lower fragment of Fig. 14 is a zoom of the respective inter- 
val of the upper fragment. We see that even such a small 
zone of the atom-field interaction is enough to generate 
a complicated self-similar structure. 

In order to prove that fractality is a robust property of 
randomly walking atoms, we compute the exit times T 
in dependence on the detuning <5 with the space between 
detectors equal to 21 wavelengths. The result is shown 
in Fig. 15a at n = 11880 photons. Comparing it with 
the same function T(S), but computed with two wave- 
lengths, we see that the structure in Fig. 15b is more 
ordered demonstrating under zooming an intermittency 
of smooth and irregular parts with singularities. 

The exit time T, corresponding to both the smooth and 
unresolved intervals, increases in average with increas- 
ing the magnification factor. There are also exist atoms 
never reaching the detectors in spite of the fact that they 
have no obvious energy restrictions to leave the small 
space between two detectors. There is an infinite num- 
ber of unstable periodic orbits and chaotic orbits. The 
measure of such orbits is, of course, zero, and they can- 
not be exactly detected even in numerical experiments. 
However, such a sensitive dependence of the measured 
characteristic on the detuning should manifest itself in 
wild oscillations of the exit times in an arbitrary small 
neighbourhood of such values of S where the scattering 
function T(S) is singular, and the range of variation of T 
should not decrease to zero as the size of such a singular 
neighbourhood is reduced. The scattering function is sin- 
gular on a Cantor set of values of the control parameter. 
It means that a relatively small uncertainty in the control 
parameters (or in initial values of dynamical variables) 
can often make prediction of the exit time practically 
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impossible. The fractal structure has been observed only 
with randomly walking atoms, i.e. at negative detunings. 
We have found complicated structures in the scattering 
function at 5 > 0, however under a magnification they 
appear to be smooth and without singularities. 

To clarify the role of dissipation in forming fractal 
structures in the numerical scattering experiments, we 
plot in Fig. 16a the same dependence T(5) but without 
dissipation, i.e. "f a = 0. The fractal in Fig. 16a is a 
Hamiltonian one, and it possesses a beautiful self-similar 
structure with intermittent smooth and singular zones. 
Comparing Figs. 16a with j a = and 15b with 7 a = 0.3, 
we see that the main part of this structure disappears 
under influence of dissipation. Instead of singularities, 
one can see smooth intervals in Fig. 15b. In order to 
demonstrate how in occurs, we plot in Fig. 16b the func- 
tion T(S) under the same conditions but with a reduced 
dissipation, 7 a = 0.2. In places of unresolved structures 
in Fig. 16a, we see some structures that, however, are re- 
solved under a magnification into smooth intervals. Wc 
may conclude that spontaneous relaxation changes the 
fractal structure of the Hamiltonian scattering function. 
It smoothes second-order self-similar structures of the 
Hamiltonian fractal Q, but their border singular points 
become fractal-like dissipative structures. 

VII. CONCLUSION 

We have performed a theoretical and numerical study 
of the phenomenon of synchronization between internal 
and external degrees of freedom of atoms in a standing 
laser wave and of related topics. The character of the 
centcr-of-mass atomic motion can be understood with 
the help of the dependence of the friction force on the 
average momentum F(|p|). Zeroes of this function arc 
either attractors or repellors. As a result, atoms tend 
to be grouped around respective values in the velocity 
space, the effect known as velocity grouping. 

Synchronization is a common rhythm of co-existence 
of mechanical oscillations of atomic center-of-mass mo- 
tion and internal Rabi oscillations that manifests itself 
in the phase space as limit cycles with different periods 
and riddled basins of attraction. Wc have derived a set of 
coupled equations for external and internal atomic vari- 



ables and obtained approximated analytical solutions for 
period-1 limit cycles as for atoms, trapped in wells of 
the optical potential, as for ballistic atoms. The fea- 
sible way to detect synchronization in real experiments 
is to measure spectra of atomic fluorescence which for a 
pcriod-m limit cycle should contain sideband frequencies, 
uif ± jkfv s /m (ballistic atoms), or UfijCl^/m (trapped 
atoms). Due to the velocity grouping, limit cycles can 
be detected in spectra of light emitted by an ensemble of 
atoms with an initial velocity distribution. 

Increasing the number of photons in the standing wave 
n, we have found cascades of bifurcations with higher 
periods to appear which eventually settle in a "stable" 
chaotic strange attractor with a positive value of the 
maximal Lyapunov exponent A and a fractional Haus- 
dorff dimension. A broadband noise, added to the equa- 
tions of motion, destroys a fine structure of the cascades, 
but the prominent period-1 and period-3 limit cycles sur- 
vive under a weak noise. 

Computing the Lyapunov map, we have found wide 
ranges of the control parameters, n and the detuning 6, 
where A > 0, and the atomic motion is expected to be 
chaotic either in the form of strange attractors in the 
wells or as flights with irregular oscillations of the veloc- 
ity. Moreover, there exists a regime where the atomic 
motion is neither one of them but resembles strongly a 
random deterministic walking with erratically jumping 
atoms. The random walking is fractal in the sense that 
a scattering function is very sensitive to small changes 
in the control parameters or initial conditions and has a 
self-similar structure with singularities on a Cantor set of 
values of one of the control parameters. We have demon- 
strated it by computing a time of exit of atoms from a 
given space of the standing wave in dependence on n or 
5. 



VIII. ACKNOWLEDGMENTS 

The work was financially supported by the Program of 
the Prezidium of the Russian Academy of Sciences, by 
the Russian Foundation for Basic Research (project 02- 
02-17796) and by the Far-Eastern Division of the Russian 
Academy of Sciences. 



[1] V.G. Minogin, V.S. Letokhov, Laser Light Pressure on 
Atoms, Gordon and Breach, New York (1987). 

[2] A. P. Kazantsev, G.I. Surdutovich, and V.P. Yakovlev, 
Mechanical Action of Light on Atoms, World Scientific, 
Singapore (1990). 

[3] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon Interaction — Basic Processes and Appli- 
cations (Wiley, New York, 1992). 

[4] S. Stenholm, Rev. Mod. Phys. 58, 699 (1986). 

[5] C.S. Adams, M. Siegel, and J. Mlynek, Phys. Rep. 240, 



143 (1994). 

[6] A. Vaglica, Phys. Rev. A 58, 3856 (2001). 

[7] VYu. Argonov and S.V. Prants, Zh. Eksp. Teor. Fiz. 

123, 946 (2003) [JETP 96, 832 (2003)]. 
[8] M.Yu. Uleysky, L.E. Kon'kov, and S.V. Prants, Comm. 

Nonlin. Sci. Numer. Simul. 8, 320 (2003). 
[9] S.V. Prants and L.E. Kon'kov, Pis'ma Zh. Eksp. Teor. 

Fiz. 73, 200 (2001) [JETP Lett. 73, 180 (2001)]. 
[10] S. V. Prants and V. Yu. Sirotkin, Phys. Rev. A 64, 

033412 (2001). 



13 



[11] S. V. Prants, M. Edelman, and G. M. Zaslavsky, Phys. 

Rev. E 66, 046222 (2002). 
[12] S.V. Prants, Pis'ma Zh. Eksp. Teor. Fiz. 75, 777 (2002) 

[JETP Letters 75, 651 (2002)]. 
[13] S.V. Prants and M.Yu. Uleysky, Phys. Lett. A 309, 357 

(2003). 

[14] A.P Kazantsev, Zh. Eksp. Teor. Fiz. 66, 784 (1974) [Sov. 

Phys. — JETP 39, 784 (1974)]. 
[15] Ch. Huygens, Horologium Oscillatorium (Apud F. 

Muguet, Parisiis, France, 1673) [English translation: The 



Pendulum Clocks (Iowa State University Press, Ames, 
1986)]. 

[16] A. Andronov and A. Witt, Zh. Prikl. Fiz. 7(4), 3 (1930). 

[17] J.C. Sommerer, E. Ott, Nature 365, 138 (1993). 

[18] S.V. Prants, in: Chaotic Dynamics and Transport in 
Classical and Quantum Systems, edited by G. Zaslavsky 
and M. Courbage (Kluwer Academic Publishers, Dor- 
drecht, NATO AST Series E, 2004). 



14 



(a) 



(b) 



400 - 




400 - 



200 - 



1 


1 !• . ■ 
















■ Ml n T 











-3.95 



-3.9 

5 



-3.85 



-3.8 



-3.95 



-3.9 

5 



-3.85 



FIG. 15: Dependence of the exit time T (in units f2 1 ) on the detuning 8 (in units of fio) with different resolutions at n = 11880 
and 7a = 0.3: (a) the space between detectors is 21 wavelengths and (b) 2 wavelengths. 



